Spectroscopy of collective excitations in interacting low-dimensional many-body 

systems using quench dynamics 
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We study the problem of rapid change of the interaction parameter (quench) in many-body 
low-dimensional system. It is shown that, measuring correlation functions after the quench the 
information about a spectrum of collective excitations in a system can be obtained. This observation 
is supported by analysis of several integrable models and we argue that it is valid for non-integrable 
models as well. Our conclusions are supplemented by performing exact numerical simulations on 
finite systems. We propose that measuring power spectrum in dynamically split ID Bose-Einsten 
condensate into two coupled condensates can be used as experimental test of our predictions. 



In many cases interesting quantum phenomena appear 
not in properties of the ground states but in the coherent 
dynamics of the system away from equilibrium. Canoni- 
cal examples from basic quantum mechanics include Rabi 
oscillations in two level systems, collapse and revival 
in the Jaynes-Cummings model, Landau-Zener tunnel- 
ing. Analysis of non-equilibrium coherent dynamics in 
a strongly correlated many-body systems is a more chal- 
lenging task than in a single particle quantum mechanics, 
hence progress has been made only in a few cases. This 
includes changing parameters in the transverse field Ising 
model [l[ , mean field dynamics in systems with BCS pair- 
ing of fermions Q, collapse and revival phenomena Q 
and non-equilibrium superfluid phase oscillations 0] of 
bosonic atoms in optical lattices. In this paper we discuss 
another example of non-equilibrium coherent dynamics 
in strongly correlated many body system. We consider 
a sudden quench of interaction parameter in nonlinear 
interacting system and show that the subsequent time 
dynamics exhibit oscillations with frequencies given by 
the poles of the scattering matrix of this many-body sys- 
tem. The sudden quench can be thus used as a tool 
for spectroscopy of the interacting models. We make 
our analysis by increasing complexity of examples: we 
start with a quantum Josephson junction (qJJ) model 
then generalize it for a Gaussian model and support the 
main observation on the the quantum sine-Gordon model 
(qSG). At the end we generalize the main statement for 
non-integrable models and make supporting arguments 
in favor of our conclusion for general models. 

The qSG model, is a popular prototypical example of 
nonlinear interacting quantum system. This model ap- 
pears as an effective description of variety of condensed- 
matter, statistical physics and field theory problems. 
Thus in condensed matter it describes low-dimensional 
spin and charge systems, disordered systems (see @ for 
review) and coupled Bose condensates @. Given the 
large freedom of tunability of parameters in ultracold 
quantum gases we have in mind the test of a theory pre- 
sented below in these dynamically decoupled ID conden- 
sates. 



The dynamics in qSG model can be analyzed using its 
exact solution. In general, describing dynamical proper- 
ties of a many-body system using the exact Bethe Ansatz 
solution is not a straightforward procedure. In our case, 
progress can be made connecting the problem of time 
evolution from a certain initial state to the equilibrium 
sine-Gordon model with a boundary. For conformally in- 
variant models a similar approach of mapping temporal 
evolution after a quench to a class of boundary phenom- 
ena was discussed in Ref. Q ■ In our case the system can 
have different mass scales corresponding to solitons and 
their bound states. 

The Hamiltonian of the qSG model is given by 



Hsg = \J dx[n 2 (x) 



(cW) 2 -4Acos(/?0)]. (1) 



In application to split condensates the <j>(x) — <f>\ ~ <f>i is 
the relative phase between the two of them and II(x) is 
the conjugate momentum proportional to the difference 
between local densities and (5 — ^2it/K Q. The inter- 
ference experiments, such as reported in Ref. Q, measure 
the 4>{x) between the two condensates. After the splitting 
the system is in a state which is not an eigenstate, hence, 
this initial state will undergo a complicated quantum dy- 
namics controlled by the many-body Hamiltonian. We 
assume that at t = Q the system is prepared in a state 
with = for all x. In reality <f> is a wave packet with 
the width determined either by the rate with which the 
condensates were separated [Io[ or by the depletion 
The Luttinger parameter K of individual condensates in 
Eq. (Q} is large for weakly interacting bosons, and ap- 
proaches 1 in the hard core repulsion (Tonks-Girardeau) 
regime. K < 1 for fcrmionic and spin systems. Generally 
K can be extracted from the microscopies 

We concentrate our analysis here on dynamical prop- 
erties of one-point correlation function, however gener- 
alization of our approach for the multi-point correlation 
functions is straightforward. The one-point correlation 
function L dx(ip(t — 0)\exp(i/3<j>(x,t)\ip(t = 0)) corre- 
sponds to the amplitude of interference fringes, A(i)Q 
and can be measured as a function of the evolution time 



2 



t. To characterize the time evolution of A(t) we analyze 
the power spectrum P(uj) = lim-r^oo | J dt e luJt A(t)\ 2 . 

It is instructive to analyze first a simpler situation, 
when the spatial fluctuations of the phase <p are energet- 
ically forbidden and the Hamiltonian |T]) reduces to the 
qJJ model 



(2) 



with E c = (3 2 hv s /L and J = 2LA where v s is a sound 
velocity. Following the standard approach for analyz- 
ing sudden perturbations in quantum mechanics we de- 
compose the initial state into eigenstates of the Hamil- 
tonian ([2]), \ip(t — 0)) = J2 n an \ n )- Eigenstates \n) 
can be found explicitly using Matthieu's functions [111]. 
Only even n's are present in decomposition of \ip(t = 0)) 
and the occupation probabilities \a n \ 2 decrease with in- 
creasing n, provided that the initial state has a finite 
width. Letting system to evolve for time t we find a 



state \4>(t)) = J2n a n e ~ 



The amplitude A(t) 



<*(t)|cos($|*(i)> = E„ m «>m(cos0) nm e- l (— »)*, 
where (cos0)„ m is the matrix element of cos0 between 
the states n and m. The power spectrum for a specific set 
of parameters is shown in Fig. [I] The inset on this graph 
shows the energy level structure and a couple of possi- 
ble oscillation frequencies. The power spectrum consists 



FIG. 1: Power spectrum for a single qJJ consisting of two 
sites. At t = the tunneling amplitude is suddenly reduced 
from a very large value to J = 10, the interaction strength 
E c = 1. The inset shows the energy levels of the qJJ with the 
arrows indicating a couple of possible oscillation frequencies. 

of different peaks corresponding to the overlaps between 
various even energy levels. The largest (central) peak in 
the power spectrum comes from beating of n = and 
n = 2 states (i.e. the (cos 0) 2O e l( " 2 ~" o)t -I- c.c. term). 
The dominance of this peak comes from a combination 
of the two factors: i) The matrix elements (cos0) mn are 
the largest for the smallest difference between n and m 
ii) The state n = has the largest weight in the de- 
composition of \tp(t = 0)). There are also contributions 
at frequencies uj m — luq for m = 4,6, ... (only m = 4 is 



shown), which give rise to peaks at higher frequencies. 
For a harmonic potential such peaks would appear at 
precise multiples of the frequency of the central peak. 
But because of the unharmonicity such peaks are shifted 
to lower frequencies. Similarly other important contribu- 
tions to the power spectrum come from beating between 
states n and n + 2 with n > 0. They appear at smaller 
frequencies than that of the central peak. Finally we 
have contributions from beatings of various other combi- 
nations of states, but they come with a smaller weight. 

We now discuss a case of full ID model ([TJ. In the 
large K limit, the qSG model can be well approximated 
by the Gaussian theory with a massive term ~ m<f> 2 , 
where m ~ J±, instead of the cosine. The initial state 
is then a product of squeezed states for all k vectors 
|^o) ~ Ilfc ex P( a ! a -fc)|0)i where a\ is the usual bosonic 
creation operator for the excitation with the momentum 
k and the energy LOk = \Jk 2 + m 2 . Expanding it one ob- 
serves that the state |?/>o) contains various combinations 
of pairs of particles with opposite momenta. The leading 
contribution to A{t) coming from two-particle excitations 
behaves at long times as (cos(</>(ir, t))) ~ sin(2mt)/Vt. 
This behavior leads to the power law singularity in the 
power spectrum: P(u>) ~ \lo — 2m\~ 1 . Higher harmonics 
corresponding to multi particle excitations correspond to 
weaker singularities and their weights are suppressed. 

With decreasing K nonlincarities start to play increas- 
ingly important role. As we will show below the peak 
at fundamental frequency ujq splits into a sharp singu- 
larity and a two-particle contribution (see Fig. [2]). The 
two-particle contribution corresponds to the excitation of 
a pair of two lowest energy breathers B\, which are di- 
rect analogues of the massive excitations in the Gaussian 
model. The singularity corresponds to the non-decaying 
excitation of an isolated E>2 breather, which is a bound 
state of two E>i breathers (for more explicit notations 
see below). In addition to splitting of the fundamental 
peak one finds that higher harmonics are shifted to lower 
frequencies in a direct analogy with a qJJ. Beatings of 
different harmonics lead to appearance of singular peaks 
at frequencies smaller than uiq, similarly to a qJJ system. 

For finite K a convenient complete basis of the 
qSG model is provided by the asymptotic scattering 
states which can be obtained by the action of the 
elements A ak (9) of the Zamolodchikov-Faddeev alge- 
bra on the vacuum state: \0\62---0 n ) ai a2 a „ = 
< (0i )A\ 2 (9 2 )...AiJ0 n )\O). Here the operators A ak (0 k ) 
have internal index a k (corresponding to solitons (+), an- 
tisolitons (-) or breathers (n) (denoted by B n below)) and 
depend on the rapidity variable 6k , defining the momen- 
tum and the energy of a single quasi-particle excitation of 
mass M a : E — M a cosh(0), P = M a sinh(0). The opera- 
tors are defined in such a way that A ak (0k) |0) = 0, Vfc, a^. 
Explicit dependence of M a (A) was found in Ref. [1 31 ] - 

Next we consider the initial condition cf>(t = 0) = 
as a boundary-in-time condition of the Dirichlet type. It 
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belongs to the class of integrable initial (boundary ^con- 
ditions, which generically have the following form [14| : 

E„ 6 ,„(^SL(0)+ 7 ff K- b (9)Al(-e)Al(9)) 

\B)=Ne - |0).(3) 

The matrix K is related to the reflection matrix of 
solitons, antisolitons and breathers. These matrices 
are known [14| for all integrable boundary conditions. 
The important feature of these boundary states is the 
presence of various poles in the reflection amplitudes 
K ab . These poles correspond to bound states of the 
particles A a ,Ab- In the soliton-antisoliton channel 
such bound states are breathers with masses Mb„ = 
2M s sia(mr/(8K - 2)), where n = 1, ...,AK - 1 and 
M s = M S (A, The poles in the soliton-antisoliton 



FIG. 2: Power spectrum for K — 1.6 and A = 0.4 including 
single and two-breather contributions. Many other contribu- 
tions are not visible on this scale. Arrows indicate <5-peaks 
of type (i) (see the text below). The inset shows the power 
spectrum from the exact diagonalization on 6 x 2 sites for the 
bosonic Hubbard model. Note the break in the vertical scale. 
Frequency units are defined by M S (A, K) for given A and K. 

reflection amplitude imply the possibility of ex- 

changing breathers with zero momentum, B a (0). The 
corresponding " H — B a " coupling is denoted by g a . For 
more discussions and results on the structure of boundary 
states see . We note that the state ([3]) is a generaliza- 
tion of squeezed states in the Gaussian theory. Presence 
of additional breather terms at zero rapidity is the con- 
sequence of unharmonicity of the underlying model. An- 
other important difference with the harmonic theory are 
the unusual commutation relations between A a 's which 
are neither bosons nor fermions. The state ([3]) is the 
initial state in our problem. The Hamiltonian (QJ is di- 
agonal in A's, which allows us to compute the time de- 
pendence of \iPb) 



\fl>{t) B ) = A/"exp 



iM b„_ t 



(4) 



a,b,n 



where M = M a + Mf,. For a translation invariant sys- 
tem A(t) = L(iP(t)\e if3 ^ x=0 ^\ip(t)} where L is a sys- 
tem size. We evaluate A(t) by expanding the exponen- 
tial form of \ij){i)B) and computing terms one by one. 
We use the form-factors (FF) technique discussed by 
us in Ref. Q for a related problem (for a general re- 
view of the FF approach see Ref. [HI)- Our calculations 
rely on recent progress in evaluation of FF for the qSG 
model [l5| . In general, FF are defined as expectation val- 
ues of the operators in the asymptotic states \9\ . . . 9 n ), 
which can always be represented in the canonical form, 
F° = (Q\Q\6± . . . 9 n ) a , . It is known from general 
principles [12J that the FF expansion converges rapidly 
with increasing number of the participating states. We 
will therefore discuss the most important contributions 
to the power spectrum, which come from the first terms 
in the expansion of Eq. (J3]). These contributions split 
into the following categories: i). Delta function peaks 
corresponding to contributions of individual breathers. 
They appear at frequencies u^ 1 ' = Mg. 2m and cl/ 11 ) = 
M B 2m — Mb 2 „ , m> n. The strengths of these peaks are 
given by g m F B . 2m /2 and g m gnFB 2m F B . 2n /4: correspond- 
ingly. Here the single-particle breather FF's are Q: 



7>exp(i^0) 



g /3 \/2sin(7rng) exp[/(-g n )]e""' 
tan (f)(cot(^)n;r i 1 cot 2 (^))V2' 



(•5) 



/OO jn 
— CTtv{-iMtcosYi{6))K ab {d)Ai{-d)AU6) 
-oo 47T 



where n = wr(l - <), £ = 1/(4K - 1), 1(9) is given 
in Q and the numerical factor Qp was computed in 
Ref. [I3|- These contributions generally decrease with 
n and decrease with K. ii) Two-particle A a Ab — contri- 
butions of excitations with equal masses corresponding to 
uj( 2 ) — 2MA a - They come from a continuum part. Their 
weights can be estimated by multiplying corresponding 
single particle FF's. iii) Interference contributions involv- 
ing more than two particles of equal or unequal masses. 

In Fig. [2] we present the results of our calculations for 
specific values of K and A. For large K the dominant 
contribution comes from beating of the vacuum state 
and the B\(ff)B-y{—6) pairs corresponding to the massive 
phonons. The corresponding contribution is absent for a 
QJJ. The reason is that B\ breathers can be excited only 
with nonzero momentum. However, the non-decaying 
peaks, corresponding to particles excited at strictly zero 
momentum, have their direct analogues (see Fig.[T]). In- 
deed the B2 — peak corresponds to the "central peak" in 
the QJJ picture. The satellite peaks B4 — B2, B§ — B4, 
B4 — etc. also have their analogues in Fig. [TJ Such 
peaks correspond to the non-decaying oscillations. Their 
S-peak nature is a consequence of the qSG integrability. 

Our analytic results are supported by exact numerical 
simulations in a system of N — 12 particles in two chains 
each of 6 sites. We used the Hubbard model with peri- 
odic boundary conditions and with the interaction U = 1 
and intrachain hopping J = 1 to realize a system with 
moderate value of K far from the Mott state . At time 
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t = the inter-chain hopping J± was abruptly decreased 
from a very large value to J± — 0.1. Despite seemingly 
small size we point that this system contains more than 
10 6 states, which is more than enough to distinguish inte- 
grable versus nonintegrable dynamics (see e.g. Ref. [HI). 
The results of such simulations are in complete qualita- 
tive agreement with the analytic predictions (see the inset 
in Fig. [2]). We can identify peaks corresponding to vari- 
ous breathers. There are also two additional two-particle 
contributions denoted as BiB\ — and B2B2 — O in the in- 
set. These peaks correspond to excitations with nonzero 
relative momentum. However unlike in the true thermo- 
dynamic limit, these oscillations are not broadened, since 
there is no continuum of different momentum states. 

Previous analysis for the integrable models can be gen- 
eralized further for non- inte grab le models using the form- 
factor perturbation theory (l9j |. If we consider not-so- 
strong deviations from integrability, which means that 
there are still no multiple particle production in a theory 
(which is equivalent to the absence of a branch-cuts in 
a physical strip of a Mandelstam variable), one can ar- 
gue that the first effect is the change of the position of 
the poles in the scattering matrix. This variation of pole 
can be shown to be given by the particle-antiparticle FF 
F(6i — 62), and the change of the particles mass is given 
by [l9| 8m 2 a ~ F® a (iir) corresponding to perturbing op- 
erator O. Therefore the leading effect of deviation from 
integrability is a shift of positions of peaks. 

It is possible to argue in general that the time evo- 
lution of expectation value of some local observable op- 
erator in initial state formed by a sudden quench will 
provides the information about the spectrum of the the- 
ory after this transition. One can indeed show [2(| that 
the boundary reflection amplitude R(k) is related to the 
bulk Green function G(x,x';t — t 1 ) as G(x,x';t — t') = 



l( e ik(w)\x-x'\ + R(k)e ik ^\ x - X '\) for the field 
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theory with dispersion k{oj). This formula is somehow 
reminiscent to the T-matrix formulation of the impurity 
scattering problem. We therefore observe that the poles 
structure of the Green function is encoded in a boundary 
reflection factor. Another part of the boundary state is 
the form-factor. The form-factors expectation values can 
be expanded using the LSZ formula[21J which relates the 
n-particle form-factors of operator O to a n-point func- 
tion. All that leads to conclusion that spectrum genera- 
tion during the quench dynamics is a generic phenomena. 

In this paper we showed that the dynamics after 
quench in many-body interacting system can be used for 
spectroscopy of collective excitations. We have explicitly 
demonstrated this on quantum sine-Gordon type models 
and argued that this conclusion is valid for more gen- 
eral, non-integrable models. We supported our statement 



for general many-body systems using connection of the 
boundary reflection amplitude to the Green's functions. 
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